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Abstract. Using N-body/hydrodynamic simulations, the influence of non-equipartition processes on the thermal 
and dynamical properties of cosmological baryonic gas is investigated. We focus on a possible departure from 
equilibrium between electrons, ions and neutral atoms in low temperature (10*-10^ K) and weakly ionized regions 
of the intergalactic medium. The simulations compute the energy exchanges between ions, neutrals and electrons, 
without assuming thermal equilibrium. They include gravitation, shock heating and cooling processes, and follow 
self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming 
coUisional ionization equilibrium. At high redshift, a significant fraction of the intergalactic medium is found to 
be warmer and weakly ionized in simulations with non-equipartition processes than in simulations in which the 
cosmological plasma is considered to be in thermodynamic equilibrium. With a semi-analytical study of the out 
of equilibrium regions we show that, during the formation of cosmic structures, departure from equilibrium in ac- 
creted plasma results from the competition between the atomic cooling processes and the elastic processes between 
heavy particles and electrons. Our numerical results are in agreement with this semi-analytical model. Therefore, 
since baryonic matter with temperatures around 10* K is a reservoir for galaxy formation, non-equipartition 
processes are expected to modify the properties of the objects formed. 

Key words, cosmology: theory - large-scale structure of the universe - intergalactic medium - galaxies: formation 
- hydrodynamics 



1. Introduction bound structures^. The hierarchical scenario causes the 

^ , IGM to separate into two phases: the hot intra-cluster 

In the basic picture of galaxy formation, baryons fall medium at T > 10^ K and the cold intergalactic medium 

into dark matter potential wells located along the net- at a few times 10"* K l)Cen et al.ll993lEvrard et alJll994j) . 

work of sheets and filaments formed by gravitational^insta- But the process of accretion and cooling of baryons inside 

Mity ("Sunvacv & Zcrdovich""l972'; |Rees fc Ostrjkcr'ig??*, bounded structures is still not well understood and a study 

piU^|1977; .White fc Rees..l97g). In this hierarchical pic- of the thermodynamic properties of the IGM is crucial to 

ture, galaxy clusters result from the large scale collapse our understanding of galaxy formation, 

of gas and from mergers of subunits. These massive struc- ^he cosmological plasma accreted in bound structures 

tures (M > 10^ Mq) are mainly present at low redshift acquires thermal energy through gravitational compres- 

(z < 2) . Mergers and accretion of matter involve hydrody- ^^^^ and shock heating. The increase in temperature heats 

namical shocks raising the intracluster medium^ tempera- preferentially heavy particles and owing to the large dif- 

ture to more than 10^ K. On the other hand in galaxy- 

size structures, shocks and gravitational compression do ^ In large scale structure formation simulations, the word 

not heat the intergalactic medium to such high temper- structure refers to any potential well created by the dark mat- 

atures and the gas can cool to a few times 10^ K and ter. This sense is then more general than referring only to 

concentrate to form galaxies. The intergalactic medium galaxies, since in such potential wells a number of galaxies 

(IGM) refers to the cosmological gas in gravitationally '^'^^ ^o™' ^ence the cosmological gas in gravitationally bound 

structures, the medium inside potential wells, is often referred 
to as intergalactic medium in the literature jCen fc Ostrikej 
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ference in mass between heavy particles and electrons the 
energy transfer between b oth species is poorly efficient 
l|Zerdovich fc Raized Il96^ . The outer regions of bound 
structures then consist of a non-equilibrium two temper- 
ature gas. In this paper we introduce non-equipartition 
processes between ions, neutrals and electrons of the cos- 
mological plasma in numerical simulations of large scale 
structure formation in order to investigate their influence 
on the thermodynamic properties of baryons. The results 
are compared to simulations in which equipartition pro- 
cesses are forced. The simulations are performed with an 
Eulerian N-Body/hydrodynamical code including gravita- 
tion, shock heating and radiative cooling processes. 

Depending on the depth of the potential wells, the 
temperature of the accreted baryonic matter ranges from 
10'' to 10* K. At high temperatures, the accreted matter 
is totally ionized behind shock fronts and departure 
from equilibrium between ions and electrons is due 
to the quite long energy exchange timescale between 
these species. The thermal decoupling between ions 
and electrons in galaxy clusters has been considered by 
several authors JFox fc Loebl |l997t IChieze et all 119981 
Takizawa fc Mineshied Il998t iTakizawal Il998|) . Indeed 
Chifeze'et'^Lrijiggsli show that temperatures of these 
two species can be significantly different, by a factor of 
~ 3, in the outer regions of galaxy clusters. Our large 
scale simulations do show such temperature differences in 
massive structures at low redshift. 

The simulations also show departure from equilibrium 
in structures much less massive than galaxy clusters. 
These structures are characterized by a warm plasma with 
heavy particle temperature between lO'' and 10^ K and 
electron temperature at a few times 10^ K. That plasma 
is also weakly ionized. We show that the fraction of warm 
IGM dramatically increases at high redshift because of 
the larger number of low mass structures, formed in hier- 
archical models. Departure from equilibrium is not only 
due to Coulomb interactions between charged particles, 
but also due to mechanisms of s hort-range forces between 
electrons and ne utral particles ijPetschek fc Bvronlll957t 
IShchekinovlll99ll) . In such low density IGM, equilibration 
timescales between these species are short compared to 
the Hubble time and we propose a semi-analytical model 
to explain them. This model suggests some insights into 
the thermodynamical story of cosmological baryonic gas. 

In this paper we focus on the low mass structure case 
and examine in detail the thermodynamic properties of 
the baryons. We show how non-equipartition processes, 
inside bound structures, at high redshift and before the 
reionization epoch, yields a large fraction of the IGM to 
be warm (between 10'*-10^ K) and we discuss the impli- 
cations for galaxy formation. The paper is organized as 
follows. Section|2]presents the numerical code and simula- 
tions. In section |31 influence of non-cquipartition processes 
on the thermodynamic properties of baryons are analyzed. 
Section ^ gives, with the aid of a semi-analytical model. 



the physical origin of departure from equilibrium in the 
low mass structure case. Finally Section |5l discusses the 
cosmological implications for galaxy formation. 

2. Simulations 

Simulations are performed with a 3 dimensional N- 
body/hydrodynamical code coupling a Particle-Mesh code 
to compute gravitational forces with an Eulerian (pseudo - 
Lagrangian) hydrodynamical code i Tevssier et alJll993) . 
The gravitational potential is solved from the dark matter 
and the gas density fields. The main feature of the code 
is the computation of non-equipartition processes between 
ions, neutrals and electrons. Each species of the cosmolog- 
ical plasma has its own specific internal energy, Ci, Cn, Go, 
satisfying the energy conservation equation: 



De 

'd^ 



-p{V ■v) + '^ -A„ 



(1) 



where px is the density, Cx the specific internal energy, 
Px = n-xkTx the partial pressure and the number den- 
sity for the X species. We do not distinguish the velocity v 
for the three species. 4' is a dissipative term due to shock 
heating and Anet is a net cooling term containing the 
radiative cooling rates (heating and cooling processes) 
and the energy exchange rates between species. We now 
briefly discuss the two terms of the energy equation, 
shock heating and elastic or inelastic processes, computed 
separately in the code. 

Shock h eating is treated using the artifi cial viscos- 
ity method l|Von Neumann &: Richtmveill95cill . To reduce 
numerical dissipation, the contribution of shocks to the 
energy is computed from a pseudo-entropy equation: 



Ds 
Dt 



(2) 



where the pseudo-entropy Sx for the species x is defined by 
exV^~^ with gamma being the ratio of specific heats and 
Vg the specific volume. The viscous dissipation function is: 

* = -gv(V-t;), (3) 

where the term equals — p(Ci£ -I- C2e^) with 

e — Ax'V ■ v/cs- Ax is the physical length of the 
cell and Cg the sound velocity. The constants Ci and 
C2 are determined experimentally and set to 1 in 
the simulations (the sensitivity to these constants are 
discussed in section 13. 3|) . The advantage of solving an 
equation for the pseudo-entropy rather than the internal 
energy is that an adiabatic flow remains strictly adiabatic. 

The net cooling term is expressed as the contribution 
of three terms: 



A, 



Ti + A + Qcxch, 



(4) 



where 7i is a heating term for example due to photoioniza- 
tion processes (see section EJ, A is the sum of the cooling 
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rates and Qcxch denotes the energy exchange between in- 
teracting species. 

Radiative coohng processes considered here are: col- 
hsional excitation, coUisional ionization, recombination, 
bremsstrahlung and Compton scattering. Ionization and 
heat input from the ultraviolet radiation background are 
not included (this point is discussed in section \^ . The 
computation of cooling rates needs to follow the chemi- 
cal evolution of a primordial composition hydrogen- helium 
plasma with nn and nue the number densities of hydrogen 
and helium, respectively. It must be noted that coUisional 
ionization equilibrium is not assumed here. The density 
evolution for the six species of the plasma (H°, He°, 11+ , 
He+, He++, e) are solved with: 



-/3hcO?Jc?IHc" + "Hc+"•cf^Ho+ 

/3Ho+f^of^Hc+ - aHo++'^c^^Hc++ 



dt ' 



(5) 
(6) 
(7) 



dt ' 

where (3i and ai are the ionization and recombination rates 
presented in Tabled and the conservation equations: 



n-^o + n-ti+ = nil, 
"-He" + "Hc+ + "-H0++ = nHc, 
nii+ + ?1ho+ + 2nfjo++ = «c- 



(8) 
(9) 
(10) 



The number densities of the six species allow computation 
of cooling rates and since the numerical code allows ions, 
neutrals and electrons to have their own temperature, 
cooling rates depend on electron temperature and are 
detailed in Table El 

Energy transfers between the three populations are ex- 
pressed by: 



(11) 



where t,j;-^^x2 is the characteristic timescale of interaction 
between xi and X2 species. The energy exchange between 
electrons and ions is d ue to Coulom b interactions and the 
exchange timescale is llSnitzeilll96i : 



To,ion(s) « 251.5 



A T, 



3/2 



no ln(A) ' 



(12) 



where A the molecular weight. To the electronic tem- 
perature and ln(A) the Coulomb logarithm. The en- 
ergy exchange between electrons and neutrals is due 
to short-range forces and, using the classical "hard 
reflecting sphere" cross-section ae,neut ~ 10~^^ cm^ 
( Drainc the exchang e timescale can be 

written as IZel'dovich fc R,aizerlll96fil) : 



'^e,neut (s) 



A 



(13) 



4.5 X 10-13 Tin Tc 

The energy exchange timescale between ions and neutrals 
is very short compared to the other timescales above. 



Table 1 . Ioni zation and rec ombi nation rate s expressed in 
s-i cm3 ffrom lBlackl lll98ll) and[Ce3 fHil) with Te„ = 

To/10" K the electron temperature. 



Ionization rates; 



/3hO = 5.85 X 10-11(1 +r^Y2)-lg(-157809.1/T„)yV2 
/^HeO = 2.38 X 10-" (1 + r,V')-^e(-2»5335.4/T.)yV2 



'g5 



X 10-12(1 + r^i/2)-ie(-63i5i5/i'=)r<!''2 



Recombination rates: 



aH+ = 8.4 X 10-" (1 + T°i')-V^^T: 
aHe+ = 1.5 X io-iOrO-6353 



aHe++ = 3.36 X 10-^° (1 + Tee 



Table 2. C ooling rates expressed in erg s ^ cm ^ (from 
iBlacd l|l98lh with modifications introduced bv lCenI l|l992h 

for temperatures exceeding 10^ K) with Ten = To/lO" K 
the electron temperature, gg the Gaunt factor taken to 
1.5, and Tbg the cosmic microwave background tempera- 
ture, Tbg ^ 2.7(1 + z). 

CoUisional excitation: 

Aho = 7.5 X 10-i«(l + r,V^)-ie(-"«3«/^«)nenHO 

AHe+ = 5.54 X 10-"(l + rjf )-^e(-4^3638/T.)y-0.397^^„^^^ 



CoUisional ionization: 



Aho = 1.27 X 10--''(l + r^ 



^In I (-157809. l/Te)^l/2 



Aneo = 9.38 X 10-^^(1 r,V')-^e(-2«5335.4/T„)yV2„^^^^^ 
AHe+ =4.95 X 10-^^(1 + r,V')-'e(-«3i"5/^"'r-o-^«VnHe+ 

Recombination: 

Ah+ = (3/2)fcTe X 8.4 X 10-''\l +T°^'')-''T^^^T^f -^71,7111+ 

AHe+ = (3/2)fcTe X 1.5 X 10-i°ri'-3'5'*^ncnHe+ 

Ahc++ = (3/2)feTe X 3.36 x 10-^°il + T°^)-^T^^%-^°-^n,nii,++ 

Bremsstrahlung: 

Aff = 1.42 X IQ-^'' gBT^^^(nii+ +71^,+ + 4nii,++)n, 

Compton interaction: 

Ac = 5.65 X 10-3«(1 + z)''(n - Tbg)ne 



It is worth noting that the term Qoxch represents an 
energy gain for electrons and an energy loss for the heavy 
particles. 

The simulations are performed for a flat A-dominated 
cold dark matter model defined by fim = 0.3, f^A = 0.7, 
0.02, h = (ffn/lOO) = . 7 The transfer func- 
tion is taken from 'Bardeen et_alJ (I198(tI ) with a shape 



t^ c 



parameter from iSueivarna. (.199^. We use th e normal- 
ization on COBE data ijBunn fc Whitelll997t) giving at 
R — 8 h~^Mpc a filtered dispersion erg = 0.91. The com- 
putational cubic volume is of 16 comoving /i-^Mpc on a 
side with iVg — 256^ grid cells and A^p = 256^ dark mat- 
ter particles allowing a spatial resolution of 62.5 h^^kpc. 
The dark matter particle mass is Mdm = 2.51 x 10^ Mq 
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Fig. 1. Isocontours of electron (left panel) and heavy particle (right panel) temperatures for a high mass structm'e 
extracted from the large scale S^t simulation including non-equipartition processes (the slices are 8 comoving /i^^Mpc 
on a side and contour levels are: 10^ 5.10^ 10^ 2.10^ 5.10"^, 10^ K). 



and the gas mass resolution is Mhm = 3.87 x 10^ Mq. In 
Section mi the simulations are 32 comoving h^^Mpc on 
a side with other parameters similar. 

We refer to simulations with non-equipartition pro- 
cesses as SsT and the results are compared with sim- 
ulations in which equipartition processes are forced. In 
that case the three species have the same temperature 
and these simulations are denoted Sn. 

3. Thermodynamic properties of baryons 

Accretion of baryons into potential wells created by dark 
matter involves gravitational compression and hydrody- 
namical shocks. Dissipative processes convert the kinetic 
energy of the cosmological plasma into thermal energy and 
the increase in temperature depends on the particle mass 
(eq. I19|l . Owing to the large difference of mass between 
electrons and heavy particles (ions and neutral atoms), 
the energy transfer between these species is poorly efh- 
cient, leading to an out of equilibrium plasma. The next 
two subsections describe the thermodynamic evolution of 
the cosmological plasma, first inside high mass structures 
comparable to galaxy clusters and, secondly inside low 
mass structures comparable to proto-galaxies. 

3.1. High mass structures or totally ionized IGM 

Figure n focuses on a high mass structure in the S3T sim- 
ulation and displays electron and ion temperature distri- 
butions. This cluster-like structure is extracted from a 
large scale field at 2: = 0. Isocontours illustrates shock 
heated baryonic matter at temperatures around 10^-10^ 
K which is typically the hot intra-cluster medium tem- 
perature. Indeed, galaxy clusters result from the merger 
of sub-units of matter inside deep potential wells created 
by the large scale collapse of dark matter. Hydrodynamic 
shocks induced by mergers heat the heavy particles to 
high temperatures and in the outer regions, just behind 



the shock fronts, the electron temperature is significantly 
lower than the ion temperature. 

In the case of discontinuities, the post-shock tempera- 
ture for each species of the plasma is computed with the 
conservation equations which connect impulsion, pu, pres- 
sure, p, and specific enthalpy, h — H/ p, on both sides of 
the thin compression region. Ions and neutral particles 
are grouped under the term heavy particles with the sub- 
script HP. Pre-shock quantities in the unperturbed region 
are denoted with the subscript and no subscript is re- 
served for post-shock quantities. The Rankine-Hugoniot 
relations, written with the compression ratio w = Po/Pi 
are: 

pu = poMo, (14) 
P = PoUo(l - ^) +Po, (15) 

h = ]^ul{l - LU^) + h^. (16) 

The specific enthalpy for the plasma is expressed by: 

/i = - ^nnpfcTHP + "^UckT^ , (17) 

with ne and nnp the number densities of the electrons 
and the heavy particles, respectively, and To and Thp the 
electron and the heavy particle temperatures, respectively. 
The increase in temperature is then calculated by assum- 
ing that the number density per mass unit for each species 
{ux/ p) is constant across the discontinuity: 

T T ^^CQ rr.^ , PoMo(l~CU^) , . 

JhP-JhPo = (Jco-Jc)H TT , (18) 

«HPo ofcnHPo 

The electron gas undergoes adiabatic compression across 
the discontinuity and th e electro n temperatu re inc reases 
by a factor of l/w^-i llZel'dovich fc Raize j Il966|) . The 
variation (Tc„ — To) can then be neglected in equation 
(|18|l . Moreover for strong shocks, 1/cj « 4, and the last 
term becomes ~ (3/16)(poI'^/(fc"-HPo))- By introducing 
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Fig. 2. Isocontours of the gas temperature in the Sn simulation (left panel) and the heavy particle temperature in the 
S3T simulation (right panel) at z — 4 (the slices are 8 comoving /i~^Mpc on a side and contour levels are: 10., 10**, 4.10'* 
K). 



the molecular weight of the gas A, the increase in temper- 
ature for heavy particles becomes: 



3 Amp 2 



2.2 X 10^ K A 



V 



100 km s" 



(19) 
(20) 



As gravitational compression and shocks heat prefer- 
entially the heavy particles, the temperature of the two 
species can differ significantly in an extended relaxation 
layer during a timescale which depends on the energy 
transfer between particles. Plasma in galaxy clusters is 
fully ionized and for typical values of density and temper- 
ature in outer regions, an estimate of the characteristic 
timescale for the energy exchange between electrons and 
ions gives (see ea. lT^ : 



9.4 X 10^ yr 



(Te/S.lO^ K)3/2 



(nc/10 



-5 



5)(ln(A)/30)- 



(21) 



This relaxation timescale is of the order of the Hubble 
time th, and explains the departure from equilibrium in 
the outer regions such as seen in Fig. ^ For higher den- 
sity, this timescale decreases and the equilibrium can be 
rapidly recovered between both species: the distributions 
in Fig. n] show that, in the structure center, both tem- 
peratures are similar. A number of earlier studies show 
such differences between temperatures in galaxy cluster 
high-resolution simulations based on Hoffman-Ribak ini- 
tial conditions ieze et al.lll998tlT^kizawa fc Mineshig3 
[l998,.: Takizawa 1998). Here we reach similar conclusions 
but obtained using simulations of large scale structure for- 
mation. 

Departure from equilibrium in the outer regions of a 
massive galaxy cluster is a first illustration of the influ- 
ence of non-equipartition processes on the IGM. However, 
in hierarchical models massive galaxy clusters exist mainly 



at low redshifts and an out of equilibrium plasma at tem- 
perature higher than 10^ K is only expected at these red- 
shifts. We now analyze the influence of non-equipartition 
processes on the cosmological plasma over the evolution of 
the universe, at epochs when a larger fraction of baryonic 
matter is contained in low mass structures. 



3.2. Low mass structures or weakly ionized IGM 

Figure |21 illustrates departure from equilibrium in struc- 
tures much less massive than galaxy clusters and at high 
redshift. Temperature distributions for a typical region 
show significant differences between the simulation with 
non-equipartition processes and the simulation in which 
equipartition is forced between species: heavy particle 
temperature in the Syr simulation is larger than plasma 
temperature in the Sn simulation. These differences are 
located in the outer regions of bound structures and are 
in the range 10^-10^ K. 

Figure 121 displays the temperature profiles along a line 
of sight through one of the structures shown on the tem- 
perature distributions in Fig.|21(the line of sight is marked 
by the dashed line). The outer parts of the profiles show 
that, in the S'st simulation, the heavy particle tempera- 
ture rises up to lO'^ K whereas the electron temperature 
is around 10'* K. In the same regions but for the S'it sim- 
ulation the plasma temperature is at a few times 10'* K. 
Then in the range of temperature 10*-10^ K, the IGM in 
the iS'sT simulation is significantly warmer than in the S'it 
simulation. 

Figure 01 plots also the ionization degree profiles along 
the same line of sight. The ionization degree is defined by 
X = nc/{n}i + Tifjo)- These profiles point out that in the 
regions of temperature differences (clearly seen on the left 
side of the profiles), the plasma is weakly ionized in the 
simulation whereas it is partially ionized in the Sn 
simulation. The former being two orders of magnitude 
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2.50 
X (h 'Mpc) 



2.50 
X (h"'Mpc) 




Fig. 3. Temperature and ionization degree profiles along 
the line of sight marked by the dashed line in Fig. |2] gas 
temperature in the Sn simulation (top panel) , heavy par- 
ticle temperature (middle panel, plus sign), electron tem- 
perature (middle panel, diamond sign), ionization degree 
in the S3T simulation (bottom panel, plus sign), ioniza- 
tion degree in the Sn simulation (bottom panel, diamond 
sign). 



lower than the latter. 

We now quantitatively estimate the mass fraction of 
that warmer plasma. First of all the baryonic mass fraction 
per interval of heavy particle temperature and per interval 
of electron temperature is computed in the S3T simulation 



at different redshifts. Figure ^ shows that the baryonic 
mass fraction with heavy particle temperature between 
10'''-2.10^ K is larger than the baryonic mass fraction with 
electron temperature in the same range. That warm out 
of equilibrium plasma exists at any redshift but its mass 
fraction decreases with redshift. 

To inquire about the density of the warm plasma, we 
plot in Fig. O isocontours of the baryonic mass fraction 
per interval of temperature and per interval of baryonic 
density contrast = Ph/{ph)- The diagram on the left, 
as a function of the heavy particle temperature, shows a 
region with temperatures in the range 2.10^-5.10^ K and 
baryonic density contrasts around 10 or less. The middle 
plot shows that electron temperature in this region is a 
few times 10^ K. The point here is that this warm plasma 
region is not found in the Sn simulation where the plasma 
is rather cold (right panel). Differences between the two 
simulations are then identified by a warm temperature 
plasma with relatively low baryonic density contrast. 

Finally we estimate the evolution with redshift of the 
different phases of the plasma in the S3T simulation and 
compare with the Sn simulation. We define a bulk plasma 
temperature in the former simulation as: 



(fie + Sj + S^ )Vs 



-2/3 



Tic + rii + Un 



(22) 



where A^tot is the total number of particles. The evolu- 
tion of the baryonic mass fraction is computed for dif- 
ferent ranges of temperature, Tg in S^t and T in Sn- 
Mass fractions are normalized to the total baryonic mass 
with a temperature higher than 9.10'^ K. Three ranges 
of temperatures are considered: 9.10'^ < T < 2.10* K, 
2.10* <T< 5.10^ K and T > 5.10^ K. Figure shows 
that, between z = 7 and z = 3, about 30-40% of the 
plasma is warm (2.10* < T < 5.10^ K) in the S3T sim- 
ulation whereas this phase represents only less than 5% 
in the Sn simulation. We note that, in particular at high 
redshift, the plasma that is warm in the Sst simulation 
is cold (9.10^ <T < 2.10* K) in the 5it simulation, and 
that the cold baryonic mass fraction is much lower in the 
simulation (this point will be addressed in the final 
section). 

3.3. Influence of the viscous dissipation function 

Before deducing any cosmological implications from the 
thermodynamic and chemical differences between the two 
simulations as described in the previous section, we have 
to check that our results are not due to numerical ef- 
fects. In fact the pseudo-entropy equation (eq.[21) computes 
the contribution to the energy of shocks and gravitational 
compression and includes a numerical heating term given 
by the artificial viscosity. Is this term responsible for de- 
parture from equilibrium seen in the S^t simulation? To 
invalidate this hypothesis, we run the same simulations 
except that now the constants Ci and C2 in the term Qv 
(eq. OJ are reduced by a factor 10. 
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Fig. 4. Baryonic mass fraction per interval of electron temperature (thin line) and per interval of heavy particle 
temperature (thick line) at different redshifts in the Sst simulation. 
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Fig. 5. Baryonic mass fraction at z = 5 per interval of baryonic density contrast and per interval of heavy particle 
temperature (left panel), per interval of electron temperature (middle panel) in the Sst simulation, per interval of gas 
temperature in the Sn simulation (right panel), light grey is high mass fractions and temperature is shown on the 
ordinate axis. 



Figure [7| plots the same temperature field as in Fig. |21 
for heavy particles in the S^t simulation and for plasma 
temperature in the S'it simulation. The temperature 
profiles are also displayed along the same line of sight as 
before. Temperature distribution in the S'it simulation 
shows that, due to the reduction of the shock contribu- 
tion, the amount of energy created by shocks is lower 
resulting in a too low pressure to support gravitational 
compression. The spatial extension of structures are then 
globally smaller than in Fig. [21 



Nonetheless, comparing the two simulations reveals 
that the S^t simulation shows also a substantial fraction 
of the IGM warmer than in the Sn simulation and re- 
gions which depart from equilibrium in the outer regions 
of bound structures are spatially larger than the region 
over which the artificial viscosity is applied. We can con- 
clude that the heating of heavy particles is not dominated 
by numerical effects. Moreover as will be pointed out in 
section l^m increasing the resolution but keeping the same 
value for the viscous term used in section leads also 
to departure from equilibrium. The thermodynamic dif- 
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Fig. 6. Evolution with rcdshift of the baryonic mass frac- 
tion computed for different ranges of temperature (nor- 
malized to the total baryonic mass with a temperature 
higher than 9.10^ K): 9.10^ < T < 2.10^^ K (plus sign), 
2.10^ < T < 5.10^ K (triangle), T > 5.10^ K (star), in the 
SsT (upper panel) and Sn (lower panel) simulations. 



a semi-analytical model. 

We consider here the relaxation of an out of equi- 
librium hydrogen plasma having just undergone gravita- 
tional compression. Relaxation processes are cooling and 
non-equipartition processes. We neglect other processes 
such as the expansion of the universe or compression due 
to the accretion of matter or merger events. The plasma 
is composed of hydrogen nuclei (neutrals and ionized par- 
ticles) and electrons with number densities nu = khp = 
nno +«H+ a-nd ric, respectively. Densities can be expressed 
with the ionization degree x : 



= xnn — rtH+ EiJ^d njjo = nii{l — x). 



(24) 



Initial conditions are suggested by the results of the simu- 
lations: temperatures of heavy particles and electrons are 
respectively, Thp = 10^^ K and = 2.10'' K, hydrogen 



nuclei density is rtu = 5.10 ^ 
degree is x = 10^^. 



cm ^ and the ionization 



We solve an equation of energy balance for heavy par- 
ticles and for electrons: 



'cxch: 



HP 



1 
--( 
P 

Qcxch ~ Qcooh 

p p 



(25) 
(26) 



where the right hand terms represent the change in the 
specific internal energy for the two fluids. Qcooi is the rate 
of energy lost by the electrons (per unit of volume and per 
unit of time): 

Qcool = Aexc(H) + Aio„(H) + Arec(H) + Abrcm(H). (27) 



ferences between the Sst and 5it simulations are then of 
physical origin and examined in the next section. 

4. Analytical interpretation 

In the outer regions of not too dense structures, the IGM 
is out of equilibrium for temperatures in the range 10''^- 
10^ K and is weakly ionized (Fig. |21 and 13) . Contrary to 
the outer regions of galaxy clusters in which the IGM is 
totally ionized and with sufficiently high density to give a 
relaxation timescale of the order of the Hubble time (see 
eg. I21|l . relaxation timescales between ions and electrons 
and between neutrals and electrons are now very short 
compared to the Hubble time : 



(23) 



We take into account interactions between neutrals and 
electrons since the plasma is weakly ionized. 

What are the mechanisms which control departure 
from equilibrium and is the relaxation timescale long 
enough to allow physical effects on cosmological time? 
This section provides answers to these questions based on 



The cooling rates associated with each process, coUi- 
sional excitation, coUisional ionization, recombination, 
bremsstrahlung, are given in Table 13 Qexch is the rate 
of energy exchange between electrons and heavy particles 
(per unit of volume and per unit of time): 



Q 



cxcli 



3, (Thp -5 

^ '^exch 



(28) 



where Tcxch includes two terms due to two interaction 
mechanisms, that between charged particles and between 
neutral and charged particles: 



1 



1 



1 



'^c,ion "^ajHout 



(29) 



We consider isochoric transformations, dp/dt = 0, then, 
since p « mpUu with TOp the atomic mass, dnn/dt = 
and equations and become: 



3 dTnp 

— knyi 

2 dt 

3, dTe 
2 dt 



Qcxch ) 

~Qcxch Qcool 2^"^6^H"^^' 



(30) 
(31) 
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Fig. 7. Influence of the viscous dissipation function: the S^t and the Sit simulations are now computed with the 
Qv term (eq. ^ reduced by 10. Upper panels: isocontours of gas temperature in the Sn simulation (left panel) and 
of heavy particle temperature in the Sst simulation (right panel) for the same field as on Fig. |21 (isocontours at 
10., 10^, 4.10* K). Bottom panels: temperature profiles along the line of sight shown on top panels by the dashed 
line: gas temperature in the 5'it simulation (left panel), heavy particle temperature (right panel, plus sign), electron 
temperature (right panel, diamond sign). 



where x is the ionization degree. The last term of the 
second equation is computed from the evolution equation 
of the electron density: 



-TT = "H-TT = Paoncnuo - aH+ricn-H+, (32) 



where /3ho and aH+ are respectively the ionization and 
recombination rates for hydrogen (see Table We 
compute the evolution of To, Thp and x with respect to 
time using equations pilj) . lpTT|l and (p^ . 

Figure shows that departure from equilibrium is 
maintained on a timescale of gigayears, and that during 
this time the ionization degree does not evolve (Fig.|S|D). 
It is interesting to note that just before the equilibration 
between the two temperatures, the electron tempera- 
ture slightly increases because of the energy gained in 



elastic processes^. Figure IHt presents the evolution of 
the different characteristic times: the energy exchange 
timescale between electrons and ions (eq.inj, the energy 
exchange timescale between electrons and neutrals (eq. 
I13|l . and the resultant Toxch due to both mechanisms (eq. 
I29|l . It is interesting to see that in the relaxation region 
and since the plasma is weakly ionized, elastic processes 
are dominated by neutral-electron interactions and 
Tcxch ~ Tc.ncut- Equilibration is recovered once plasma 
is getting partially ionized and during this short pe- 
riod, Toxch ~ To. ion- Oncc thc ionization degree increases, 
Toxch decreases drastically and equipartition is established. 

Moreover, Fig.|Sli shows that the energy exchange rate 
between heavy particles and electrons is approximately 
equal to the cooling rate. This means that there is no 
change in the electronic internal energy (eq. 126(1 . By writ- 

^ Si milar evolutions we r e com puted in a sophisticated anal- 
ysis bv iFadevev fc Gilleli il998l) . of the structure of radiative 
shock waves propagating in stellar envelopes. 
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(c), and energy exchange and cooling rates (d) using equations (120)), and 



ing Qexch ~ Qcooi we then derive a relation between Thp 
and Te such as: 

THP=re + (|^)-lQcool. (33) 

Since Qcooi = /(^o, ?t.o, njjo, this relation is only pa- 
rameterized by the ionization degree'^. In Fig.lOt- we plot 
Thp evaluated in equation (|33|l . with the cooling term 
given by equation (|27|l . for different values of the ioniza- 
tion degree. In a narrow electronic temperature range, 
10^-3.10^ K, the heavy particle temperature extends in 
a much larger range, between lO" and 5.10^ K. It IS in- 
teresting to note that the higher the value of the ion- 
ization degree, the smaller is the difference between the 
two temperatures. For very low ionization degree the re- 
lation between Tjjp and To does not depend on x since 
7o,nout oc {n-naTl^'^)^^ , and Fig. IHt shows that departure 
from equilibrium in weakly ionized plasma is controlled by 
neutral-electron mechanisms. 

In this range of electronic temperatures, cooling via 
bremsstrahlung processes is not efficient, and neither is 
cooling due to recombination processes (too low density). 
Then only terms due to coUisional excitation and coUi- 
sional ionization processes are effective in equation (|27|) . 
Figure illustrates the relation between Thp and 
when equation is solved by neglecting the coUisional 

^ This approach was used by IPetschek &: BvronI (^S^) to 
determine the ionization rate of a gas of argon in a shock tube. 



excitation term. We note that even in that case, departure 
from equilibrium is also seen, but now for higher electron 
temperatures. 

Then for a given ionization degree and electron 
temperature slightly higher than 10"* K, the energy 
exchange rate between the heavy particles and the 
electrons compensates for the energy lost by the electron 
gas due to inelastic processes. We now return to the S^t 
simulation and confirm this scenario by plotting temper- 
ature values for heavy particles and electrons in each cell 
of the computational volume. Each symbol in Fig. [TUl 
corresponds to the range of ionization degree of each cell 
of the simulation. We overplot relation (|33|l for different 
ionization degrees. For instance squares indicate cells with 
ionization degree between 10~^ < x < 10~^, and we see 
that these cells are in the region delimited by the curves 
X = 10^^ and x = 10^^. The agreement for each range of 
ionization validates the condition Qcxch ~ Qcooi- These 
results confirm our understanding of the physical origin of 
the non-equilibrium IGM in moderately dense structures, 
although numerical simulations include a larger set of 
physical processes (as expansion and accretion). 

Departure from equilibrium results from the competi- 
tion between elastic processes in the weakly ionized IGM 
and cooling processes. It is clear that the thermodynamic 
history of accreted gas in bound structures is not the 
same in a simulation with non-equipartition processes 
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Fig. 9. Heavy particle temperature versus electron tem- 
perature for different values of the ionization degree com- 
puted by equation (|33|l with (a) and without (b) collisional 
excitation cooling term in eq. H27|l . 



and in a simulation in which equipartition is forced. As 
the plasma undergoes gravitational compression in not 
too dense structures, heavy particles are heated much 
more than the electrons. The electron temperature does 
not increase sufficiently to partially ionize the medium 
(see also the evolution of temperature and ionization 
degree in Fig.lHJ. Therefore the ionization degree profiles 
in the outer regions in the Szt simulation (Fig. I^J show 
a weakly ionized plasma with Tc ^ 10^ K. But in the 
Syy simulation the electrons can be heated up to 2.10^ 
K, and though the ionization temperature of hydrogen is, 
with / the ionization potential, T = I /k ^ 1.6 x 10^ K, 
a pla sma at T 2.10^ K is ~ 90% ionized ijShchekinovl 
Il99lh . The ionization degree of a hydrogen plasma in 
collisional ionization equilibrium is /3ho/(/3h'' + ctH+)- 
This fact explains the two orders of magnitude difference 
between the ionization degree profiles plotted for the Syr 
and the S'it simulations (at 1.5 x 10^ K, the ionization 
rate for hydrogen is ~ 10"^'^ s~^ cm'^ whereas at 2.5 x 10^ 
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Fig. 10. Heavy particle temperature versus electron tem- 
perature at z = 3 in 300 cells randomly chosen in the Syr 
simulation; each symbol represents the ionization degree 
in each cell and the curves are computed in the same way 
as in Fig. 



K the ionization rate is 



cm^). 



Before closing this section, we would like to mention 
an interesting point. In the center of structures such as in 
Fig. [3 the IGM is cold and partially ionized in the S'it 
simulation (Fig. This suggests that in such not too 
dense structures the IGM can be never heated to the so- 
called virialized temperature and can be partially ionized 
and radiate its energy just by being heated up to 3.10'* 
K. This is supporting evide nce to the i mport ant ques- 
tion recently summarized by iKatz et aD l)2002j) concern- 
ing results from different numerical simulations of large 
scale structure formation, showing that a non-negligible 
fraction of the ac creted gas can never reach the virial- 
ized temperature ()Katz fc GunnI Il99lt iKav et all l2000t 
iFardal et al1l20flll) . 

5. Cosmological implications 

The introduction of non-equipartition processes alters, for 
not too dense structures, the properties of plasma with 
temperature in the range 10^-10^ K. This plasma is ex- 
pected to contribute to galaxy formation. In this section 
we qualitatively analyze the influence of non-equipartition 
processes on galaxy formation. We discuss when and how 
this influence is dominant. 



5.1. At high redshift 

The analytical interpretation shows that the weaker the 
plasma is ionized the larger the departure from equilib- 
rium. Figure El shows that the difference between heavy 
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Fig. 11. Baryonic mass fraction at z = 10 per interval of the baryonic density contrast and per interval of the bulk 
temperature in the S^t simulation (left panel) and per interval of gas temperature in the Sn simulation (right panel) , 
both simulations including photoionization processes. 
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Fig. 12. Collapsing redshift and mass of (1-2-3)(t fluctua- 
tions (from left to right), computed by a spherical top-hat 
collapse model (eq. in a A-cold dark matter universe. 



particle and electron temperature (in the range 10^-10^) 
is larger for lower ionization degree and Fig. |S1 shows that 
decoupling is mainly due to interaction mechanisms be- 
tween neutrals and electrons. This suggests that two con- 
ditions are required for a not too dense plasma to be out 
of equilibrium. The first one is that the plasma accreted 
into structures must be weakly ionized. The second one is 
that cooling processes must involve atomic cooling (coUi- 
sional excitation and/or coUisional ionization). This sug- 
gests that the major influence of the non-equipartition 
processes occurs before the end of the reionization epoch, 
i.e. at epochs when the cosmological plasma is not fully 
ionized and when the cooling is dominated by atomic cool- 
ing. 

To numerically confirm this point, we introduce pho- 
toionization processes in the S^t and the Sn simulations 
to see what happens after complete reionization of the 
universe. As the analysis of the influence of photoioniza- 



tion processes is not the aim of this paper, we do not give 
the details of the numerical aspects of this implementation 
(see Alimi & Courty, in preparation). Photo- ionization 
rates and heating rates are computed from the evolution 
of the hydrogen and helium densities and from the spec- 
trum of the ult r aviole t radiation background. Following 
IWeinberg et alJ l)l997j) the photoionizing background is 
assumed to have a redshift-dependent evolution, with 
the bulk of the transition between z = 7 and z — 6 and 
the peak between z = 3 a nd z = 1 . More over, as shown 
in the numerical study bv ICnedinl l|20n(]l) . we take into 
account a shallow evolution of the ultraviolet background 
intensity before z ~ 7 and begin the ionization processes 
at the redshift of 11.5. The radiation is considered 
as a spatially uniform field over the computational 
volume. We compute the baryonic mass fraction per 
interval of baryonic density contrast and per interval of 
temperature, which is the bulk temperature for the S^t 
simulation (eq. I22|l . Figure [TTI shows that at z = 10 a 
fraction of the plasma is warmer in the simulation with 
non-equipartition processes than in the Sn simulation. 
At z = 8 no differences can be seen between the two 
simulations. This suggests that the IGM is warmer in a 
simulation with non-equipartition processes only before 
the universe has been completely reionized. However to 
draw any firm conclusions concerning the temperature 
differences between heavy particles and electrons before 
this epoch, simulations including a more careful treatment 
of how reionization proceeds should be performed, which 
is currently out of reach. 

As mentioned above, non-equipartition requires struc- 
tures such that cooling is dominated by atomic processes 
and the question is now to examine whether such struc- 
tures exist at high redshift. Figure [T^ illustrates the mass 
and the collapsing redshift of a v-a fiuctuation computed 
in the top-hat collapse spherical model and with cosmo- 
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Fig. 13. Heavy particle (left panel) and electron (right panel) temperature distributions at z — 12 for the high 
resolution Sst simulation with Np — Ng — 384'^ (the slices are 2.166 comoving /i~^Mpc on a side and contour levels 
are: 2., 3., 10., 10^, 1.1 x lO'^, 1.7 x lO'^, 2.5 x 10*, 2.10^ K). 
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Fig. 14. Temperature profiles along a line of sight (dashed 
line in Fig.ll3|l in the structure at z ~ 12 for the high res- 
olution S^T simulation: heavy particle temperature (plus 
sign), electron temperature (diamond sign). 



logical parameters of our A cold dark matter model. Here, 
V is the relative amplitude of fluctuation of scale M in 
units of the filtered dispersion cr(M, z = 0) at z = 0: 



rcoll 
"lin 



cr(M,z = Q)D{zy 



(34) 







with the critical spherical ovcrdcnsity at z 

and D{z) the growth factor (Padmanabhan,.199^ . Before 
z = 10, which is likely before the end of the epoch of 
reionization, 3cr fluctuations have already collapsed. These 
curves are plotted in a simplified way and we refer to 
iHutchings et alJ l)2002(l for recent studies. In a virialized 
approach, structures with mass higher than a minimal 
mass of around 10* Mq are known to cool efficiently by 
atomic processes and the so-called virialization tempera- 
ture is larger than 10*^ K. But as numerical simulations 
do not make any assumptions about the temperature of 
the IGM, we perform a high resolution simulation taking 
non-equipartition processes into account to check if struc- 
tures with an out of equilibrium plasma are present at 




Fig. 15. Isocontours in a electron temperature-ionization 
degree diagram of cooling timescales computed with equa- 
tions (P5|l (dashed curves) and (solid curves) normal- 
ized to the dynamical timescale fea. lTTjl . the labels on each 
curves are the ratios icooi/^dyn- 



redshift higher than the epoch of reionization. The num- 
ber of dark matter particles and grid cells is now equal 
to 384"^ for the same computational box length, allowing 
the initial density fluctuation spectrum to include smaller 
scale modes. 

This high resolution simulation is then able to catch 
bound structures at very high redshift. An example of such 
a structure at z = 12 is illustrated in Fig. 1131 and temper- 
ature profiles along an horizontal line of sight through the 
structure is plotted in Fig. 1141 Similar differences between 
the electron and heavy particle temperatures can be seen 
in the outer regions of the structure as previously reported 
in the lower resolution simulation. It has to be mentioned 
that the computational size is L = 16 Mpc and that a 
lower box size would show a more resolved structure. We 
conclude that the introduction of non-cquipartition pro- 
cesses leads to an out of equilibrium IGM, and that their 
influence is likely to be dominant at epochs before com- 
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plete reionization of the universe. This plasma is involved 
in the galaxy formation process and we now turn to the 
cosmological implications. 



5.2. On galaxy formation 

Galaxy formation is closely related to the ability of gas to 
cool. One necessary condition for the collapse of gas clouds 
is that the cooling timescale must be less than the dynami- 
cal timescale. Cooling rates depend on the electron density 
and we have seen that, in out of equilibrium regions, the 
electron density can be lower by up to two orders of mag- 
nitude than in the simulation with equipartition processes 
forced. We compare the cooling timescales computed for 
an out of equilibrium plasma (subscript NE) and for a 
plasma in equilibrium (subscript E). In the first case the 
plasma temperature is the bulk temperature, mainly the 
heavy particle one: 



the dynamical timescale: 



(^cooi)nE 
(^cooi)e 



3/2fcntotT'g 

Qcooi(-^C: '^) 

3/2kntotT 

Qcool 



(35) 
(36) 



Note that in the second expression, the electron temper- 
ature in the cooling rate term is also the temperature of 
the plasma, T = T^. These timescales are normalized to 



t. 



dyn 



37r 

32Gp' 



(37) 



and are plotted in Fig. 1151 as a function of the ionization 
degree and the electron temperature, since this tem- 
perature controls the cooling rates. In these expressions 
Qcooi is computed from equation H27|l . ritot is fixed 
at 5.10~^ cm~'^, nn is given by ntot/(l + x) and p is 
nYirrip + neme. In equation (|37|l p is multiplied by a factor 
of 10 to include the dark matter component. 

Two points matching the thermodynamical state of 
the plasma in the two simulations are worth mentioning. 
A point located at x = 5.10""* and Tc = 1.5 x lO'' K 
mimics the warm weakly ionized baryonic matter in the 
outer regions of not too dense structures (such as in 
Fig. 121) seen in the Sst simulation. A point located at 
X = 10~^ and Tc = 2.10'* K mimics the cold partially 
ionized baryonic matter in similar structures in the S'it 
simulation. Figure 1151 shows that the first point is in 
the region icooi > ^dyn whereas the second point is in 
the region icooi < ^dyn- This suggests that the out of 
equilibrium plasma in the S^t simulation is already 
cooled in the S'it simulation. Thus the cooling timescale 
is expected to be longer in simulations taking into-account 
non-equipartition processes. 
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Wc therefore estimate the cold baryonic mass fraction 
in the two simulations. Figure 1161 presents the baryonic 
mass fractions computed at two redshifts as a function of 
the baryonic density contrast in three ranges of tempera- 
ture: 8.10^ < T < 1.7 X lO'' K, 1.7 X 10'' < T < 3.10^ 
K and 3.10"* < T < 5.10^ K. First of all we notice 
that the plots for the S^t simulation show a warm phase 
(3.10"* < T < 5.10^ K) which is not found in plots for 
the Sit simulation. As explained above, this phase in the 
S3T simulation is cold in the latter simulation and this is 
the second point to notice: the mass fraction of cold gas 
(8.10^ < r < 1.7 X 10^ K) is lower in the Sst simulation 
than in the Sn simulation. We have previously shown the 
evolution in redshift of the warm and cold fractions in Fig. 
ini Since this cold phase is a reservoir for galaxy forma- 
tion, this analysis of the thermodynamic properties and 
distributions of baryonic matter, strongly suggests that 
modifications in galaxy formation are to be expected in 
models with non-equipartition processes. This issue will 
be quantified in a forthcoming paper focusing on simu- 
lations including non-equipartition processes and galaxy 
formation (Alimi & Courty, in preparation). 

6. Conclusions 

We have considered in this paper the influence of non- 
equipartition processes on the thermodynamic properties 
of baryons in the universe, using numerical simulations 
of large scale structure formation. In the S^t simulation 
each species of the plasma, electrons, ions and neutral par- 
ticles, has its own kinetic energy. The results have been 
compared with the Sn simulation in which equipartition 
processes are forced. In the S^t simulation, plasma is out 
of equilibrium in the outer regions of structures where 
the accreted baryonic matter is shocked or gravitation- 
ally compressed. The causes of the departure from equi- 
librium depend on the interaction mechanisms and thus 
on the thermodynamics of the accreted plasma. In the 
intra cluster medium of massive structures, baryons are 
shock heated to very high temperatures (~ 10^ K) and 
the equipartition timescale between ions and electrons is 
of the order of the Hubble time. But in low mass structures 
we have pointed out that departure from equilibrium in 
weakly ionized accreted plasma is driven by interactions 
between electrons and neutral particles. In that case the 
temperature differences are in the range 10'*-10^ K and in 
the S3T simulation the plasma is found to be warmer than 
in the S'it simulation. We have checked that the thermo- 
dynamic differences are not due to artificial heating. Since 
the energy exchange timescale between neutral particles 
and electrons is relatively short compared to the Hubble 
time, we have estimated semi-analytically the relaxation 
timescale in weakly ionized low temperature out of equi- 
librium regions. Wc have shown that departure from equi- 
librium results from the competition between the cooling 
rates of the electrons and their heating rates by heavy par- 
ticles. The relaxation timescale has been found to be long 
enough to yield departure from equilibrium. Moreover a 



semi-analytical relation between heavy particle and elec- 
tron temperatures parameterized by the ionization degree 
has been derived and our numerical results are in a good 
agreement with this relation. Finally, we have discussed 
the cosmological implications of the warmer intergalactic 
medium in simulations taking non-equipartition processes 
into account: we have concluded that galaxy formation is 
expected to be modified and this issue will be studied in 
detail in forthcoming papers, in terms of galaxy properties 
and galaxy clustering. 
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